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Within the model of a Lorentz violating extension of the Maxwell sector of the standard model, 
modified light propagation leads to a change of the resonance frequency of an electromagnetic cavity, 
allowing cavity tests of Lorentz violation. However, the frequency is also affected by a material- 
dependent length change of the cavity due to a modified Coulomb potential arising from the same 
Lorentz violation as well. We derive the frequency change of the cavity taking both into account. 
The new effects derived are negligible for present experiments, but will be more pronounced in future 
tests using novel resonator materials. 
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I. INTRODUCTION 

Einstein's special relativity and the underlying prin- 
ciple of Lorentz invariance are among the foundations 
of modern physics. In the past, they have been tested 
frequently, with increasing precision. Among others, op- 
tical means (e.g., [l|,0,y]) are used for such tests. Today, 
space-borne instrumentation (SUMO Q and OPTIS p, 
for example) as well as terrestrial experiments 0,11 HI us- 
ing ultrahigh precision optical and microwave techniques 
are performed or proposed to explore the validity of this 
fundamental theory even further. One reason to con- 
tinue testing SR is because it is one of the pillars of mod- 
ern physics. Another reason is that most approaches to- 
wards a quantum theory of gravity such as string theory 
and loop gravity predict Lorentz violation at some level 

nam mill. 

Electromagnetic radiation has provided the first 
glimpse of Lorentz invariance in the famous Michelson- 
Morley experiment 0, 0] , establishing the direction in- 
variance of the speed of light c. Kennedy-Thorndike tests 

0. establish the invariance of c under observer Lorentz 
boosts. These experiments have been performed with in- 
creasing accuracy, today using electromagnetic cavities 
instead of interferometers. They are based on the modi- 
fication of c arising from Lorentz violation (depending on 
the observer frame, direction of propagation, and polar- 
ization), which in turn changes the resonance frequency 
v = nc/(2L) of a cavity of the Fabry-Perot type, where 
n = 1,2,3,... and L is the cavity length. (This applies 
for cavities where the radiation propagates in vacuum, 

1. e., the index of refraction is 1.) A change of v can be 
detected sensitively due to the potentially very high fre- 
quency stability 8v/v ~ 1CU 15 possible with cavities. 



Tests of Lorentz invariance can be analyzed within an 
extension of the standard model developed by Kostelecky 
and co-workers 0, 0] . The classical test theories, such 
as the Robertson-Mansouri-Sexl test theory [13, ll2l LU| 
or the c 2 formalism 0] , can be recovered as special cases 
of this model 0]. The standard model extension starts 
from a Lagrangian formulation of the standard model, 
adding all possible Lorentz violating terms that can be 
formed from the known particles together with Lorentz 
violation parameters that form Lorentz tensors. In the 
purely electromagnetic sector, these parameters are given 
by a tensor (kF) K .\u.v which has 19 independent compo- 
nents [TH |20) . These Lorentz violation parameters can 
be viewed as remnants of Planck-scale physics of an un- 
derlying fundamental theory, such as string theory, at 
an attainable energy scale. Cavity tests are interesting, 
since they can in principle access all 19 components of 
(kF)K\fiv- The relative change of v is calculated in 
assuming L =const, so this can be interpreted as a modi- 
fication of t he phase velocity of light in vacuum c = Cq+Sc 
with El 
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[(N x E*)(k H b)w>(N xE) — E*{n DE ) uh E] . (1) 
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Here, N is a unit vector pointing along the length of 
the cavity and E is a unit vector vector perpendicular to 
N that specifies the polarization (Fig. QJ. The asterisk 
denotes complex conjugation, (khb)^ and (nDE)\nb are 
3x3 matrices resulting from a decomposition of (kp) K Xfiv 
in the laboratory frame ^J. Comparison of Eq. (JJ to 
the experiments leads to upper limits on some of the 
(^f)kA^ in the range of 10 -8 - -15 . Because these values 
are small, here and throughout the paper it is sufficient 
to work to lowest order in the components of {kp)^^. 

However, a modification of the Coulomb potential 
<5>(x) arising from Lorentz violation in general affects the 
length L of the cavity and thus the cavity resonance fre- 
quency. This has already been pointed out in the litera- 
ture EElIlllilll. Since 5v/v = 5c/c-5L/L, a length 



2 



change SL connected to $(x) might alter or even com- 
pensate the Lorentz- violating signal of the experiments. 
Owing to the different electromagnetic composition of 
different solids, SL will depend on the material. Indeed, 
very early Michelson-Morley experiments have been re- 
peated using exotic materials (pine wood and sandstone) 
for the interferometer to exclude that an accidental can- 
cellation was responsible for the once puzzling null result 

M 

Within the standard model extension, modifications of 
c as well as $(a?) = e 2 /(47r|x|) + V(x) can be treated on 
a common basis. For many-particle systems, like solids, 
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[l5j . Here, e is the electron's charge and e a an integer, so 
that ee Q is the charge of the ath particle. The vector x a b 
is the displacement between the ath and the 6th charge. 
The factor 1/2 corrects the double-counting of pairs, the 
prime denotes that a — b is excluded from the summa- 
tion. This equation shows that in the case of a modified c 
as given by Eq. QJ, modifications of the Coulomb poten- 
tial $(x) necessarily arise. This is due to the covariancc 
of the Maxwell equation, i.e., a modified wave equation 
implies a modified Poisson equation. 

In this work, the influence of this modified Coulomb 
potential to the length of the cavity is derived. As a 
result, we find a somewhat increased sensitivity of ex- 
periments to some of the parameters (kF) K \ fMl ,. That in- 
crease is not more than a few percent for materials com- 
monly used (quartz and sapphire), but more pronounced 
for ionic materials. As a main conclusion, we thus con- 
firm that it is possible to neglect effects due to Lorentz 
violation in the Coulomb potential in the interpretation 
of present cavity tests of SR. 

In Sec. Ill Al we derive the length change SL for ionic 
crystals, giving explicit numbers for cubic crystals in Sec. 
Ill Til The case of covalent bonds is discussed in Sec. Ill CI 
The result is discussed in lllll In the Appendix, we show 
that in our model, the repulsive potential between the 
atoms in the lattice is almost not changed due to Lorentz 
violation, so we need to consider only the modification of 
the Coulomb potential. 



II. LENGTH CHANGE OF CRYSTALS 
A. General considerations 

The hypothetical Lorentz violating length change of 
a cavity depends on the material it is made of. As a 
first model, we will treat ionic crystals; we discuss the 
generalization of the result to covalent bonds later. 

A lattice is formed so that equilibrium between the 
attractive and repulsive forces between the atoms is 
reached. In ionic crystals, attraction is caused by 
Coulomb interactions between the ions. The repulsive 



potential is due to the overlap of their orbitals. Both the 
attractive and the repulsive potential may be changed 
due to Lorentz violation. In the appendix, we will find 
that the modification of the repulsive potential is negli- 
gible compared to the modification of the attractive po- 
tential, so that it is possible to assume non modified re- 
pulsive forces for the rest of this paper. 

For this analysis we consider a cube, which is small 
compared to the dimensions of the solid under consider- 
ation, but large compared to the lattice, so that boundary 
terms can be neglected. Without Lorentz violation, the 
cube's side-length is L; Lorentz violation changes this to 
L + SL 1 (the superscript i — 1,2,3 denotes the spatial 
components in the laboratory frame) , thus distorting the 
cube. The solid consists of a number of these cubes, and 
the fractional change SL % jL of the solid's length are the 
same as the fractional change for a single cube. It is thus 
sufficient to consider a cube of side-length L. The total 
energy of the lattice depending on SL 1 can be written as 



E(L + SU) = const + -EyL(5L 1 ) 2 + ^-SL l 
' 2 v ' dL l 



(3) 



where Ey is the Young modulus |25| . The term propor- 
tional to Ey is Hooke's law (for simplicity, we restrict 
ourselves to isotropic elasticity). Whereas L minimizes 
E(L) without Lorentz violation, L + SL 1 is obtained by 
minimizing the total energy including V as given by Eq. 
(J2J). In that equation, the summation ^ afc can be re- 
placed by summing over a only and multiplying with the 
total number of interacting charges N, i.e. the number of 
ions in an ionic lattice. Denoting (x l ) a b the i-ih spatial 
component of x a b, and the vector components of the 

primitive translations lj, we have (x % ) a b — n J lj with a 
vector n— (n 1 ,^ 2 ,^ 3 ) £ R 3 ; thus 
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The summation is carried out over all ions. For comput- 
ing the length change, we need the derivative 
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We have s nl = s ln and t klno is invariant under permuta- 
tions of all indices. 
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The length L of the sample is a linear combination 
of the primitive translations l 1 -. Thus, the derivatives 
dl^/dLi can be obtained. We express 
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and insert this into Eq. 10 • From that, 8L{ can be 
obtained by minimizing the energy E given by Eq. 
We will do that explicitly for a cubic crystal in the next 
section. 



B. Cubic symmetry 

To give an explicit result we consider a crystal that 
without Lorentz violation has cubic symmetry (NaCl 
structure). Lorentz violation will in general break this 
symmetry. For the crystal without Lorentz violation, 
= a5j with a lattice constant a. Since next neighbors 
have opposite charges, en — (_i) nl +™ 2 +« 5 we nnc i 
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The summation is carried out over {Z 3 \0}. t abcd is only 
nonzero if the indices are two pairs of equal indices. Thus, 

j.aabb ^.abab j.abba ^ab ^j^J^ 



t ab = 



a? {—m 



(„l)"l+"2+"3»2 n 2 
(n k n k ) 5 / 2 



T±~e , a^b, 



(-1) 
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Numerical summation (summing from —50... 50) gives 
a w -0.58, tj. w 0.23, and r|| = -1.04. Using these 
equations, we obtain (no summation over n) 
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Let the z axis be parallel to N [as it was defined below 
Eq. JU]. Thus, we have 



SU _ dV dll 
dlZ ~ dlfdL, 



8V1_ 

81% rj z 
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where the (constant) rj z = L/a is the number of unit 
cells making up L with the lattice parameter a. Setting 
to zero the derivative of Eq. we obtain the fractional 
length change in z direction. Substituting 



N = N r . 



N A Vp 

' M 



where Na is Avogadro's constant, V the crystal's volume, 
p its density, M its molecular weight, and the factor N m 
is the number of atoms (ions, in an ionic crystal) per 
molecule, we obtain 
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and 
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(kde)zz, (kde)± = (kde)xx + (KDE)yy 
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Here, we included a dimensionless factor i 2 , which mea- 
sures the effective charge of the ions. For most crystals, 
l = 1, for non ionic materials, < i < 1 (see Sec. Ill Cll . 
We also included factors v\ and v%, the number of valence 
charges for the atoms. Using the dreibein N,E,N x E, 
this gives 



5L 

— = a\\N(K DE )i ah N 



(16) 



+a ± [E*( KDE ) lah E +{Nx E*){ KDE ) lah (N x E)] . 

Here, dp = A(2a — 3r||) and a± = — 3Atj_. This is now a 
vector equation which will hold in any coordinates (this 
does not mean that the cavity length depends on the 
polarization). For this equation, it is not necessary to 
make some assumption about the orientation of the cav- 
ity axis relative to the lattice, since cubical lattices are 
isotropic. For lattices of lower symmetry, the parameters 
in Eq. fTHjl would, however, become dependent on the 
angle between the crystal axes and the cavity orientation 
N. 

In practical experiments, the material might not be 
a single crystal, but a non crystalline material. These 
consist of a large number of microscopic crystals that 
are randomly oriented. If these microscopic crystals are 
cubical, however, Eq. I|ltj|) wm hold for any of them 
as well as for the macroscopic body, since in this case 
the equation makes no reference to the orientation of the 
lattice. 

The coefficients a\\,a±, and A for some materials are 
given in Table [I] Without Lorentz violation, the ionic 
model predicts the correct distance of ions within a few 
percent error. Thus, we can expect a similar precision 
for the Lorentz- violating length change of ionic materials. 
At present, such materials are uncommon for cavities in 
high precision physics, but their use has been proposed 
|26j because of the availability of ultra-pure specimens. 

Materials of other than cubic structure will show 
a length change of similar magnitude: The structure 
enters this computation via the constants a and r. 
These are in analogy to the Madelung-constants a of 
solid-state physics, which are of order unity and do not 
depend very much on the actual crystal structure (see 
e.g. j2?J for some values). This also holds for a and 
t. However, lattices with lower symmetry will lead to 
a more complicated expression than Eq. i|14[l . i n which 



(13) 

also non diagonal terms of (k-de)^ enter. 



Length change for microwave resonators 

Sometimes, especially in microwave experiments, cav- 
ities are used where the electromagnetic radiation of the 
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cavity mode goes along the circumference of an annulus 
or a sphere. For cubic crystals, it is straightforward to 
generalize our result to a cavity of any geometry. There- 
fore, we parametrize the closed path of the mode by 
means of a parameter s. We denote the length of the 
resonator mode path L m . Its relative change is given by 
integrating Eq. i|16|) along the closed path of the mode: 



SL 
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r L a\\N(nDE)ubN + a±E*(K DE )\ ah E 

+a ± (NxE*){ KDE ) lah (NxE)]ds (17) 

with S = <f ds. Here, the dreibein N,E, and N x E 
is used, where N gives the local orientation of the wave 
vector at each value of s and E _L N gives the polariza- 
tion. This equation holds for cubic lattices (that have no 
preferred axis); it also holds for uniaxial crystals (such as 
sapphire, for example), as long as the mode is orthogonal 
to the crystal axis for all values of the parameter s. In 
this case, however, the values of aj_ and ay have to be 
computed specifically for the uniaxial crystal. 

For example, for a cavity where the mode has the ge- 
ometry of a ring with the z axis being the symmetry axis, 
we identify the parameter s as the angle < </> < 2tt (Fig. 
Pi; we have N = (cos 0, sin0, 0). Choosing E = (0,0, 1) 
parallel to the z direction, 

= 11 "t ~ [(K DE )xX + (K DE )yy]+a±(K DE ) ZZ (18) 

for the path length of a microwave cavity mode like it is 
shown in Fig. ^ This result is valid for a cavity made 
from a cubic crystal; for a uniaxial crystal, it is valid as 
long as the crystal c axis is parallel to the z axis (provided 
that a± and ay are computed for the particular material). 

In whispering gallery resonators, a dependency of the 
index of refraction on Lorentz violation would also have 
to be taken into account. Compared to the length change, 
this probably is a small correction, since in whispering 
gallery resonators only a part of the electromagnetic field 
energy travels within the material. 



C. Covalent bonds 

The materials quartz (used, for example, in 0,0) and 
sapphire (e.g., 0, IS EH) are based on covalent bonds. 
These bonds have a partial ionic character, i.e., the ef- 
fective charges of the atoms are Lev (as before, v denotes 
the number of valence charges), where i < 1 can be de- 
termined roughly from the difference of the electroneg- 
ativities of the atoms (see, e.g., 29]). From Eq. l|T5fl. 
we thus obtain length-change coefficients in the per cent 
range (Table P) . Since the concept of partial ionic char- 
acter is not an exact concept, these values have limited 
accuracy. They suggest, however, that the relative length 
change SL/L of a cavity made from quartz or sapphire 
is negligible compared to the relative change Sc/c of the 
speed of light. 



N 




FIG. 1: Left: Optical (Fabry-Perot) cavity. The unit vector 
N is parallel to the cavity axis, the unit vector E is the po- 
larization of the electromagnetic radiation inside the cavity. 
Right: Microwave cavity. The ring denotes the electromag- 
netic cavity mode, parametrized by an angle < <f> < 2n. The 
polarization E is parallel to the z axis, the unit vector N is a 
tangent denoting the direction of the mode for each cj>. 



Here, we give some arguments indicating that the 
length change of a covalent crystal may be derived in 
analogy to ionic crystals. Covalent bonds are given by 
Coulomb interactions between delocalized electrons, and 
between the electrons and the atom cores. For our cal- 
culation, the main difference to ionic bonds is that the 
electrons are not localized; instead, they are described by 
a symmetrized n-electron state ip n {x i, X2, ■ ■ • x n ). This 
state is formed from the single-electron states ipi as the 



Slater determinant (1/V7V!) Det(i/ , i(xj))- The covalent 



potential Vc = V E c + V EE 
core-core interactions Vcc 



Vcc is then the sum of the 



(19) 



[where a and b enumerate the cores, and ^ ab \x a — Xb) is 
the Coulomb potential between the cores a and b at the 
locations x a and Xb], the electron-core interactions 

fe=^^f (rr<s> {ta) (x t - x a )^ n d 3 Xl . . . d 3 x n 

i a J 

r , (20) 

(where i numerates the electrons and &^ a > is the 
Coulomb potential between the core a and the electron 
i) and the electron-electron interactions 

Vee = J2 [ & n )*® m (xi ~ x ] W l <Px l . . . d 3 x n 

(21) 

with $w, the Coulomb potential between the electrons 
i and j. Since tjj n is symmetrized, these expressions con- 
tain all the terms usually found in the treatment of cova- 
lent systems, like the hydrogen molecule. Although the 
electron wave-functions are spread over the whole crys- 
tal volume, their centers of mass are localized at some 
point within a distance a from the cores, where a is of 
atomic dimensions |30j , smaller than the typical distance 
of atoms in a lattice. As a toy model, we thus assume 
centers of charge at points x a at some distance a from 
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TABLE I: Length change coefficients. The first part of the 
table contains cubical two-atom crystals, for which the calcu- 
lation was performed. The length-change coefficients A, a«, 
and a± are obtained from Eqs. 11411 and 11511 . The second 
part contains crystals of other structure with non ionic bind- 
ings. Quartz and sapphire have two lattice constants; the axis 
of highest symmetry is given. 



Material 


E Y 


P 


M 


l?V\V2 


a 


A 


a | 


a ± 




(GPa) 


(cm*) 






(A) 








NaCl 


40 


2.1 


58 


0.64 


2.841 


-0.14 


-0.28 


0.10 


LiF 


64.8 


2.64 


26 


1 


2.015 


-0.54 


-1.06 


0.37 


NaF 


64.8 


2.73 


32 


1 


2.135 


-0.43 


-0.83 


0.29 


sapphire 


497 


4.0 


102 


1.5 


13.0 


-0.02 


-0.03 


0.01 


quartz 


107 


2.2 


60 


0.72 


5.4 


-0.06 


-0.11 


0.04 



the locations of the nearest cores. Coulomb forces be- 
tween these centers of charge and the cores make up an 
important part of the covalent binding force. The model 
reduces Vec and Vee to a sum over Coulomb interac- 
tions, i.e., we have 

V EC ~ ^^]$ (aa) (f a -x a ) (22) 

a a 

Vee ~ ££*< a5 >(af B -5r 5 ). ^ 

a j, 

The distance a between the centers of mass of the de- 
localized electrons and the cores gives rise to an effec- 
tive charge of the cores. A calculation of the Lorentz- 
violating length change could now start from these sums 
and proceed as for ionic crystals. While this is a sim- 
plified model, it indicates that the effect of a Lorentz 
violating Coulomb potential for covalent bonds will be 
in analogy to ionic crystals. Thus, estimates for covalent 
crystals made above should at least give the correct order 
of magnitude. 

III. CONCLUSION AND SUMMARY 

The contributions from Eq. (JTJ and Eq. I|16|l give 

8v - a 

— = -a«N{K DE ) uh N 
v " 

' ' ' r E*( KDE ) lah E (24) 

+ (N x E*)( KDE ) llib {N x E) . 
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This has been simplified by noting that astrophysical 
tests lead to (khb) = —{k-de) with an accuracy many 
orders of magnitude higher than attainable in cavity ex- 
periments [l5| . The second part of Eq. I|24l) is equiva- 
lent to the speed of light change Eq. (Q times 1 + 2a ±, 
a± ~ 0.01 . . .0.4 for the materials in tabled Addition- 
ally, the frequency is sensitive to a term proportional 



to an ~ —0.03 ... — 1. For both coefficients, the small 
coefficients in the per-cent range are for covalent ma- 
terials (quartz and sapphire); although they are based 
on a simplified model, they should be accurate enough 
to conclude that the relative length change of a cav- 
ity made from quartz or sapphire is negligible against 
the change of the speed of light. Since present experi- 
ments use these materials, it is indeed possible to ana- 
lyze them as if the cavity length would not be affected 
by Lorentz violation. This holds for cavities of any struc- 
ture, including (microwave) cavities using radial or whis- 
pering gallery modes. Their length change is given by 
Eq. H17fl . However, future experiments using resonators 
out of ionic materials, as proposed in [23 will be affected 
more strongly by the Lorentz violating Coulomb forces 
(see Table HJ; here, it is necessary to analyze them using 
Eq. (OJ rather than Eq. JTJ • 

We note that a cavity experiment is mainly sensitive 
to the parity odd coefficients [15| of the standard model 
extension, since both the change of the phase velocity 
of light, Eq. QJ, as well as the cavity length change 
Eq. (|16fl are dominated by these coefficients. Parity 
even coefficients that are not tightly constrained by 
astrophysical birefringence measurements |15| enter 
cavity experiments only suppressed by /3i a b, the velocity 
of the laboratory with respect to the frame in which one 
defines the parameters kde and khb- For experiments 
on Earth and using a sun centered frame, /3i a b ~ 10 . 

We have calculated the length change of electro- 
magnetic cavities due to a modified Coulomb potential 
arising necessarily from a Lorentz violating velocity of 
light, as it follows from an extension of the standard 
model. As a first model, we have assumed that the 
cavity is made of an ionic crystal. We then extended the 
model to the covalent bonds of practical materials using 
the approximate concept of partial charges. Taking 
this into account for the interpretation of cavity tests 
of Lorentz invariance, we have shown that the length 
change effect adds to the hypothetical Lorentz-violating 
signal derived from a change of the speed of light 
and increases somewhat the sensitivity of experiments. 
However, for materials used in present experiments, 
the increase is at the percent level and thus negligible, 
as all the Lorentz invariance tests are null tests. As 
a result, the classical interpretation that neglects the 
length change is justified (this also holds for microwave 
cavities and for the classic interferometer tests). For 
future experiments using ionic crystals as a resonator 
material, however, the effects derived here are larger and 
must be taken into account. 

As an outlook, we note that the properties of mat- 
ter and thus the dimensions of cavities also depend on 
Lorentz-violating terms in the fermionic equations of mo- 
tion, i.e., a modified Dirac equation 0, |^|. Strong 
upper limits on some of these terms have been placed 
by, e.g., comparisons between atomic clocks |32j. In this 
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work, we have focused on the purely electromagnetic sec- 
tor; the effect of the fermionic terms in cavity experi- 
ments will be treated elsewhere. It might be possible to 
access Lorentz violating fermionic terms in cavity exper- 
iments due to the length change they cause. One would 
compare cavities made from different materials, in order 
to separate the electromagnetic from the fermionic terms. 
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APPENDIX: THE REPULSIVE FORCE OF THE 
IONIC BOND 

Here, we calculate the modification of the repulsive 
force due to Lorentz violation. It is caused by the modi- 
fied wave functions of the ions arising from the modified 
Coulomb potential between the nucleus and the outer 
electrons. The repulsive force is a short-range force, so it 
is sufficient to consider next neighbors only. 

Since in ionic crystals, the positive ions are usually 
small compared to the negative ones [2j| , we assume as a 
model that the positive ions are pointlike. The repulsive 
potential 



$bm(R) =const|V (0) (i?)| 2 



(25) 



is then proportional to the probability density of the 
outer electrons of the negative ions. ■0 (R) is the unper- 
turbed wave function at the location R of the positive ion 
without Lorentz violation (we assume that the ion-ion in- 
teraction does not change appreciably the wave function, 
which should be satisfied to reasonable accuracy in ionic 
crystals). The wave function tp^(R) is proportional to 
e~CI R l for large \R\ with a constant £ ~ 1/ao with the 
Bohr radius ao; the potential $bm is thus a Born-Mayer 
potential £?e~ 2 ^ r (B as well as £ are phenomenological 
parameters that are determined by fitting the model to 
the measurements), which is frequently used as a model 
of the repulsive force j2?J • In the case of Lorentz viola- 
tion, ip(r) — tp^ + with a correction ipw that is 
proportional to kde, which will be calculated for hydro- 
genlike atoms. 

For this calculation, we denote the unperturbed state 
of hydrogen \nlm >(°) with a principal quantum number 
n, angular momentum I and magnetic quantum num- 
ber to. The corresponding wave function Vw«,(^D can 
be factorized into a radial function R n i(r) that depends 
solely on the radius coordinate r, and spherical harmon- 
ic) . Only the ground state n — 1 is not degen- 



ics Y, 



erate; for the n = 2 states, the use of the perturbation 
theory of nondegenerate states is still possible since, due 
to the nature of the matrix elements calculated below, 
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. 001 
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-0 .002 
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FIG. 2: Ion-ion potentials for n = 1,1 = 0, m = (for other 
quantum numbers, a similar situation is found). The x axis 
gives R in units of ao, the y axis the potential in arbitrary 
units. Solid line: 3005"1>sm(^), i.e. the correction to the 
Born-Mayer Potential, assuming {k,de) zz = 1, enlarged by a 
factor of 300. Dots: V(R), i.e. the correction to the Coulomb 
interaction between two ions for comparison. Dashes: Un- 
perturbed ion-ion Potential <I>bjv/ — 1 / (4tt 1 7?| ) . The constant 
const= 15 has been inserted as to obtain a minimum of the to- 
tal unperturbed ion-ion potential around R=5, as commonly 
found in ionic crystals (e.g., R = 5.3ao for NaCl). The correc- 
tion to the repulsive potential is negligible compared to the 
correction of the Coulomb potential. 



the degenerate n — 2 states do not mix. For n = 1,2, 



the state \nlm >= \nlrri > 



(o) 



ilm >( 1 ) in hrst order 



perturbation theory is given by 

' < nlm\V\n'l'm' > 

n'l'm' 



\nlm> {1) = J2 



E„ — E„> 



\n'l'm! > (26) 



with V given by Eq. J2J). Inserted into Eq. lt2*5|) . this 
leads to a correction to the Born-Mayer potential 



5$bm{R) = 2 const V" 



>< nlm\V\n'l'm' > 

E„, — E„> 



(27) 



n'l'm' 

*\^L(r)€Lar)\ 



We apply coordinates in which kd e is diagonal with the 
diagonal elements K xx ,K yy , and n zz . Applying polar co- 
ordinates with the 9 = axis parallel to the z axis, the 
Lorentz violating correction to the Coulomb potential V 
can be written as 

e 2 

V = [(k xx sin 2 <p + K yy cos 2 <p) sin 2 9 + k zz cos 2 9] . 

Since in first order perturbation theory, the change of 
the wave function \nlm >^ for a sum of perturbations 
Vi is a linear combination of the changes for each 
single perturbation Vi, it is sufficient to consider the term 
proportional to k zz . In configuration space, 



<nlm\r 1 cos 2 9\n'l'm' >= / R n i(r)-R n , v (r)r 1 dr 
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x / Y l m {9, </>) cos 2 9Y{? (9, 0) sin 9d9dcj> (29) 



with the normalized radial functions 

where Lz (z) are associated Laguerre polynomials as de- 
fined in |33||. For the treatment of hydrogenlike atoms, 
£ = Z/aQ is given by the core charge number Z. Here, 
we take C — l/a . Using 




and zero otherwise. That means, the matrix element is 
nonzero only for I = I' or I — I' ± 2. However, for / = /', 
the radial integral is nonzero only for n = n'; that term is 
not needed for this calculation. Combining these results 
leads to the more simple expression for the correction of 
the Born-Mayer potential 



5<S>bm(R) = 2 const V" 



< nlm\V\n'l + 2m> 
E„ — E n i 



n'— 3 



(35) 



Y? 



! 2l + l(l-m)\ 
47r (l + m)\ 



[where P[ n (x) are the associated Legendre functions], it 
is easy to verify that the (f> integration gives ir5 m ,m'- For 
the 9 integration, we substitute x = cos 9. Applying the 
relation [^3] 

x{2l + l)PT(x) = (l-m+ l)P^(x) + (/ + m)P?Mx) 

(32) 

two times and using the normalization 



it can be shown that 



4i 2 -l+(i+m)(4i 2 +4m-3) 2(l+m)\ 
(2i-l)(2Z+l) 2 (2i+3) (l-m)\ 

2 (Z+m+l±l)! 

(2Z'-l)(2Z'+l)(2Z'+3) (l-m- 



P! n (cos9)e m4 ' , (31) (for n = 1,2, < I < n, and |m| < I). Formally, 

the series should be summed up to infinity. However, 
the terms decay rapidly with increasing n' . The correc- 
tion of the Born-Mayer potential due to Lorentz violation 
can thus be obtained and compared to the correction of 
the Coulomb potential due to the same Lorentz violation 
(Fig. [21). Here, R was assumed parallel to the z direction. 
If R is not perpendicular to the z direction, the correction 
to the Born-Mayer potential is smaller. From this cal- 
culation, it can be concluded that the Lorentz-violating 
modification of the repulsive potential 63>bm is negligible 
compared to the modification of the Coulomb potential, 

(33) at least for the n — 1 and 2 states that have been con- 
sidered in this simple model. This is due to the small 
magnitude of the matrix elements < nlm\V\n'l'm' >, 
which are also small for high n and n' (they even de- 
crease with increasing n and n'). So, although for n > 3, 
the perturbation theory for degenerate states should be 
applied, this will not dramatically change the magnitude 
of 6$bm and thus our conclusion that 5$bm is negligible 

(34) remains valid. 



Pr(x)Pir{x)x 2 dx 



V = I 



L±L > ! I' — 1 + 2 

1±1)! ' ' 31 ^ 
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